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ABSTRACT: A shear deformation theory including residual thermal and moisture 
effects is developed for the analysis of either symmetric or unsymmetric laminates with 
mid-plane edge delamination under torsion loading. The theory is based on an assumed 
displacement field which includes shear deformation. The governing equations and 
boundary conditions are obtained from the principle of virtual work. The analysis of the 
[90/(±45) n /(+45) n /90] s ECT mode III test lay-up indicates that there are no hygrothermal 
effects on the mode III strain energy release rate because the laminate, and both 
sublaminates above and below the delamination, are symmetric lay-ups. A further 
parametric study reveals that some other lay-ups can have negligible hygrothermal 
effects even when the sublaminates above and below the delamination are not 
symmetric about their own mid-planes. However, these lay-ups may suffer from 
distortion after the curing process. Another interesting set of lay-ups investigated is a 
class of antisymmetric laminates with [±(0/(0-90) 2 /0)] n lay-ups. It is observed that when 
n takes on even numbers (2 and 4), both hygrothermal and mode I effects can be 
neglected. From this point of view, these lay-ups provides a way to determine the mode 
III toughness between two dissimilar layers. However, when n takes on odd numbers (1 
and 3), both hygrothermal and mode I effects may be strong in these lay-ups. In 
particular, when 0 equals 45°, the lay-ups are free from both hygrothermal and mode I 
effects irrespective of n. 

KEY WORDS: laminated composites, fracture toughness, delamination, strain energy 
release rate, mode III fracture toughness test, torsion. 

Nomenclature 

a delamination length, m 

Aij extensional stiffness coefficients, N/m 

b specimen width, m 



2?ij coupling stiffness coefficients, N 

Dy bending and twisting stiffness coefficients, Nm 
Eh Young's moduli, GPa 

Fx applied tension force, N 

g 2 strain energy release rate parameter due to mechanical effect, Nm 

g i strain energy release rate parameter due to hygrothermal and 

mechanical coupling effect, N 

go strain energy release rate parameter due to hygrothermal effect, N/m 
G T total strain energy release rate, N/m 
h half the laminate thickness, m 

AH percentage moisture weight gain, weight % 

I integration constant, m 

Mb applied bending moment, Nm 

Mt applied twisting moment, Nm 

{M} vector of stress couples, N 

n number of repeats of a group of layers 
{N} vector of stress resultants, N/m 

Qj 3-D transformed reduced stiffness, N/m 2 
s characteristic root, 1/m 

AT temperature change from cure temperature to test temperature, °C 
x, y, z Cartesian coordinates 

Greek Letters (Symbols) 

{a}, { 72} lamina coefficients of thermal expansion in lamina and laminate 
coordinate systems, respectively, per °C 
{PJ, { P } lamina coefficients of moisture expansion in lamina and laminate 
coordinate systems, respectively, per % by weight 
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£o laminate mid-plane extension strain 

£ij strain tensor 

{£} vector of in-plane strain components 

<|)i extension-twist coupling, m 

<|>2 bending-twist coupling 

{y} vector of out-plane shear strain components 
r\ normalized constant term in Af xy , 1/m 

cp l =cp 2 delaminated portion of laminate width, or a , m (see Eq. A-28) 
tp 3 undelaminated laminate width, or b-a> m (see Eq. A-28) 

X laminate torsional stiffness, Nm 2 

Ai strain energy parameter due to hygrothermal effect, Nm 

Ao strain energy parameter due to hygrothermal effect, N 

k bending curvature, 1/m 

fi ij shear moduli, GPa 

II total strain energy per unit laminate length, N 
0 twisting angle per unit length, or twist, 1/m 

0 HT twist due to hygrothermal effects, 1/m 

0 fiber angle with respect to x axis for a ply in the laminate 
Cij stress tensor 

{a} vector of in-plane stress components 

{x} vector of out-plane shear stress components 

['Ey] coefficients that relate the end tension, bending and twisting 
moments to extension strain, bending curvature and twist Nm 2 
{ x Pj} HT force and moments due to hygrothermal effects, Nm 2 
Ei, E 2 extension-twist and bending-twist coupling indicators, respectively 
¥** Vy rotation functions in warping and transverse directions, respectively 
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Introduction 

Delamination is a common failure mode in laminated composites, and often the 
most serious. Characterization of a laminated composite material for its resistance against 
such failure has become an important task for damage tolerance design purposes. 
Delamination in composites is characterized in the fracture mechanics sense by the three 
modes: mode I (opening), mode II (in-plane shear) and mode III (anti-plane shear). In 
general, the stress field at the delamination front could be a mixture of all three loading 
modes. Characterization of fracture resistance, or fracture toughness, in each individual 
mode for laminated composites can help a designer in materials selection and as a 
foundation for fracture mechanics criterion used to design against delamination failures. 

Analyses and test techniques for mode I delamination fracture toughness have 
been established [1], and resulted in the publication of an ASTM standard test method D 
5528-94a [2]. Mode II fracture toughness characterization [3-7] has also received 
extensive attention, and standardization development is under way. In the past, the 
study of mode III fracture toughness had been largely ignored due to the fact that the 
mode I fracture toughness is much lower than mode III in the traditional laminated 
composites with epoxy matrices. The lack of a suitable mode III test method at the time 
was also a factor [8]. As the mode I fracture toughness improves with the development 
of enhanced resin systems, the two shear failure modes may become more important. 
Furthermore all the three failure modes are important in the establishment of a mixed 
mode fracture criterion for structural design purposes. 

Recently, an Edge Crack Torsion (ECT) test was proposed as a mode III 
toughness test method [9]. This test is based on a class of [90/(±45) n /(+45) n /90] s lay-up 
specimen with a mid-plane free edge delamination subjected to torsion. The ECT test 
was analyzed based on a shear deformation theory, for its validity as a mode III test [10]. 
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In addition, a simplified solution for the [90/(±45) n /(+45) n /90] s ECT lay-ups was 
obtained, and a parametric study on the shear modulus effect was also performed [11]. 
The analytical solution developed in Ref. 10 was restricted to symmetric laminates. 
Asymmetry in sublaminates (lay-ups above and below the delamination) was considered, 
but the residual thermal and moisture effects were neglected. The present work is based 
on the shear deformation theory developed in Ref. 10, but includes the analysis of 
asymmetrical laminates. The hygrothermal effects are also included. The present solution 
should be identical to the solution of Ref. 10 on laminate and sublaminate symmetric lay- 
ups. The present solution will be able to identify the hygrothermal effects neglected in 
Ref. 10 when the sublaminates are not symmetric. In addition, the present analysis can 
deal with asymmetric lay-ups as well. In the following, a shear deformation theory 
including thermal and moisture effects is developed for the analysis of arbitrary 
composite lay-up with mid-plane edge delamination under torsion loading. The theory is 
applied to the ECT lay-ups first. Next, a study of some altered ECT lay-ups is presented 
to assess the thermal and moisture effects on these lay-ups. Finally, the theory is applied 
to a class of antisymmetric lay-ups. 

Mathematical Model 

A symmetric or asymmetric laminate with mid-plane free edge delamination is 
subjected to torsion loading as shown in Fig. 1 . The delamination is assumed to run 
through the entire length of the laminate. The laminate is modeled as three sublaminates 
as shown in Fig. 2. Sublaminate 1 and 2 represent the portions of the laminate above and 
below the delamination, while sublaminate 3 represents the undelaminated portion. 

D i s pla ce ment Fi eld 

The laminate is assumed to be loaded at its ends in such a way that the stresses do 
not vary along the x axis. The resultant of the stresses at a typical cross section reduce to 
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a twisting moment as shown in Fig. 1. If the laminate is loaded by a static equivalent 
twisting moment at the ends, the stresses do not depend on the x axis at locations away 
from both ends by virtue of the St.-Venant principle [12]. Such a state of deformation fits 
the definition of the so called generalized plane deformation [13]. By following the 
derivation of Ref. 13, the displacement field may be written as 


u(x, y, z) = £ 0 • x + K • x • (z + z) + K x • x • (y + y ) + U 0 (y, z) 
v(x,y,z) = -jKi-x 2 + G-x-(z + z) + V 0 (y,z) 
v/(x,y,z) = -jfC-x 2 -Q-x-(y + y) + W 0 (y,z) 


( 1 ) 


where u, v, and w denote displacements along the x, y, and z axes, respectively. The 
axial extension is denoted by e 0 , while the bending curvatures in the x-z plane and x-y 
plane are denoted by K and K lt respectively. The angle of rotation per unit length about 
the x-axis, or twist, is denoted by ©. 


The constant z represents the z-coordinate distance of sublaminate mid-plane 
with respect to the delamination plane; z=0 for sublaminate 3 because the sublaminate 
mid-plane coincide with the delamination plane, z=h/2 for sublaminate 1, and z--h!2 for 
sublaminate 2. The constant y is the distance between the twisting center and the 
delamination tip and hence is identical for all three sublaminates. The twisting center is 
initially unknown, and y does not show up in the expressions of torsional stiffness and 
total strain energy release rate of the laminate as will be seen in the final analysis. 


The bending curvature k x is neglected because the thickness-to-width ratio of the 
laminate is much smaller than unity. Further simplification is achieved by considering 
first order terms only in the Taylor expansion of the displacement functions U 0 and V 0 
with respective to the reference surface, and neglecting thickness (z) dependence in 
function W 0 as the following: 
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dU 0 (y,0 ) , 2 s ] 

U 0 (y,z) = U 0 (y,0)+ ¥■ z + Olz 2 )+-- 

0 Z 

d V 0 (y,0) / 2 \ 

V 0 (j,z) = V 0 (y,0) + — z4-0(z 2 ) + -.- (2) 

w , 0 (y.d = w r o(y.o)+^)+- 

Neglecting the terms of 0(z 2 ) and higher for U Q , V 0 and O(z) and higher for W 0 in 
Equation 2 and substituting the rest into Equation 1 yields 

u (x,y,z) = £ 0 'X+K-x-(z + z) + U(y) + z- V x (y) 

v(x,y,z) = 0-x*(z + z) + V(y) + z- yr y (y) > (3) 

w(x,y,z) = -|K-jc 2 -0-:r-(y + y) + W(y) 

Shear deformation is recognized through the rotations \| r x and \| f y . Displacement 
functions U, V , W, y* and xfy are functions of y only. These displacement functions will 
be determined from the governing equations and boundary conditions in section 
following the next section on constitutive relationships. 


The strains corresponding to this displacement field can be written as 



where 
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e° xx =£ 0 + KZ K a = K r°y Z =¥y + Wj 

£ yy = V .y K yy=Vy,y T& = Wx -©-(y+JO ■ 

fxy=U y +Q-Z ^=^+0 


( 6 ) 


where partial differentiation by a variable is denoted by the subscript comma followed 
by that variable. 


Constitutive Relationships 

A generic sublaminate along with stress resultants and moment couples is shown 
in Fig. 3. The stress resultants and moment couples are denoted by N x , N yj N^, Q x , Q r 
M x , M y and M ^ as shown in Fig. 3. The constitutive relationships can be written in terms 
of stress resultants and moment couples and associated strains and curvatures as follows 
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( 14 ) 


(W HT , {M} HT ) = |^ 2 [e]({«}Ar + {^}Atf) (l, z)-dz 
Q- = the transformed reduced stiffness [14] 

{a}= (a xx a yy a xy j l lamina coefficient of thermal expansion in x-y 
coordinates 

{/3}= (p xx ~p yy ~p xy ) T lamina coefficient of hygroscopic expansion in x-y 
coordinates 

AT- temperature differential from cure temperature to test temperature, and 
AH- percentage moisture weight gain. 

For sublaminate 3, the integrations in Equations 9, 10, 13 and 14 should be carried out 
from -h to h. 

Governing Equations 

Because all the strains and stresses are independent about the x axis, a unit length 
of the laminate along the x axis is considered. For all three sublaminates, the upper and 
lower surfaces are stress free. The governing equations and boundary conditions for the 
three sublaminates can be obtained from the principle of virtue work: 

Jj n i 4 H dydz+ IL CT u 5 4 d y dz + JL» °« s 4 dydz = F > 5u ' ds (15) 

where O 1 , Q 2 , Q? represent the volumes of sublaminate 1, 2 and 3, respectively, while T 1 , 
T 2 and T 3 represent their respective boundaries. In order to distinguish similar quantities 
of the sublaminates, superscripts 1, 2 and 3 are used to refer to the respective 
sublaminate. In the event of a quantity raised to a power, the quantity will be bracketed. 

By following the procedures presented in Ref. 10, the governing equations 
obtained from Equation 15 for each of the sublaminate may be simplified as 
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( 16 ) 

(17) 


Afj=Nj,= M\ =C;=0, j =1,2,3 
M iy, y -Qi= j = 1.2,3 

The boundary conditions and continuity conditions are 


Ml A = 0 

y \y=-{b-a ) 

(18) 

Ml A = J = 0 

1 y—a 1 y=a 

(19) 

m ’,L- m :,L + <L 

(20) 


(21) 


Solutions for the Displacement Functions 

Substituting the stress resultants and moment couples given in Equations 7 and 8 
into Equation 16, and using Equation 6 yields 


.^ + 0z4 = -[« j ]j Co l-{^}-(< y + 0)+{« j }, j = l,2,3 (22) 

. <y I 

v 4 +w .i = -7r[^ -©-( t+t)]. j =1.2,3 (23) 

— . A , 

where [R } ], {/? J } and {fl J } are given in Equations (A-l) through (A-3) of the Appendix. 


The remaining non-zero stress resultants and moment couples can be simplified as 
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A ^ 

The matrix [§>] and vectors { Q] and{ $} are given in Equations (A-4) through (A-6) of 
the Appendix. 


The solution of \| f x for each sublaminate can be obtained by applying Equation 17 
to the three regions respectively and making use of the boundary conditions given in 
Equations 18 through 21. The general form of can be written as 


Y i x = I i -(e siy + e 2s ' p ' ■e- s,y )+t-e s ‘ p ‘ ■e- s ‘ y + @-(y+y), j = l,2,3 


(26) 


where V represents the integration constant, p ] is replaced by -(b-a) for sublaminate 3, 
and by a for sublaminate 1 and 2, respectively. The integration constant is different from 
each sublaminate. The three constants for the three sublaminates are given in Equation 
(A-7) of the Appendix. 


The characteristic root s ] and parameter in Equation 26 are given by 


s ] = 


j _ 


'Ms 


■^44 



j =1,2,3 


t 

n j = iii e o+iilKT+20+2-, j = 1,2,3 


(27) 

(28) 


Loading Conditions 

The extension force, bending and twisting moments for the entire laminate cross- 
section are denoted by F x , A/ B and A/ T , respectively. These are expressed as 


Fx = j.V a) ^-^+J 0 a K + ^ 2 )-^ (29) 

M b = Ml ■ dy + \^M\ + Ml +z'-N\+z 2 - N 2 X }dy (30) 
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( 31 ) 


m t = 2 C (b - a) M i- d y + 2 l 0 l M iy +M w + z l -K +z 2 ' N l\ d y 

Substituting the stress resultants and moment couples from Equations 24 for the three 
sublaminates into Equations 29 through 3 1 yields 


F x 

M a 

J?L 

n 

* & 
-V ' 

Mt. 


M 



UJ = 1,2,3] 


(32) 


The parameters 'Ey in Equation 32 are functions of sublaminate stiffness coefficients (Ay, 
5y, Dy) and specimen geometry ( a , b , h ) and are given in Equation (A-18) of the 
Appendix. The parameters 'F J Hr represent the thermal and moisture contributions, and 


their expressions are given in Equation (A- 19) of the Appendix. 


The strain e 0 , bending curvature k and twist 0 can be written as a combination of 
the mechanical and hygrothermal contributions as 
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► +< 
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where 


HT 


K 

0 


HT 


r i-l r -i HT 

=-N W 


(33) 


(34) 


For a laminate under torsion loading, F X =M B = 0, the induced extension strain and 
bending curvature can be written as 

e 0 =<p v Q+<^ r 

K=<j>2 •© + 02 ^ 
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(35) 

(36) 



where 



'F 


1-1 


12 


22 


¥ 
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13 


23 


''Pa 

*12 ' 

"1 

fv.l 

L 21 

^22. 

1 

l*aj 


HT 


(37) 

(38) 


Parameters <{>! and <t> 2 represent extension-twist and bending-twist couplings, 
respectively. 


The twisting moment can be obtained from Equation 32 as 


M r = A-(©-0 HT ) (39) 

where 

A — (f)y • O 3 1 + (f>2 ' ^32 ‘^‘^33 (40) 

The coefficient X represents the laminate torsional stiffness. 


Strain E n ergy a nd S trai n Ene rg y Re l ea se R a te 

The strain energy per unit length of the laminate can be written as 

n = ^X(Lj ({£}-{«} A 7 ’-{^} Afl') T {cr} dydz + \ Ayf {?}dydz) ( 41 ) 

1 j=l ' 

Substituting Equation 11 into Equation 41, rearranging yields 

3 3 

n = ^X Lj({e} T {cr} + { r } T {T}) d'y • rfz - X Lj <{«}-Ar+{^} A^) T [G]{f} dydz 
1 j=l 1 j=l 

3 

+ \t,L{m^T+{p}AH) T [Q]{{a}AT+{P}AH) -dydz 
1 j=l 

(42) 

The last term in Equation 42 does not depend on delamination length a , hence will not 
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be considered further. The first term in Equation 42, could be easily proved to be 


n,=-(e 0 ■F x + kM b +QM t ) (43) 

For a laminate under torsion loading. Equation 43 reduces to 

n, =i©M T =—A©(0-e HT ) (44) 

^ 2 

The second term in Equation 42, n 2 , may be written as 

n 2 = -I©A 1 -i.A 0 (45) 

Parameters A x and are functions of delamination length, a, and given in Equations (A- 
26) and (A- 27) of the Appendix. 

For the laminate (as show in Fig. 1) containing a planar delamination of length a 
that extends under a constant twist 0, no external work is performed as the 
delamination extends, the total strain energy release rate may be calculated from 

G t =-£2- (46) 

o a 

Substituting Equation 44 and 45 into Equation 46 and differentiating yields 

G T = -^h a ® 2 +^(h a Q m + + A 1>fl )*0 + -A Oa (47) 

In order to isolate the hygrothermal effects, the total strain energy release rate can 
be written in terms of mechanical twist by substituting Equation 33 into 47 to yield 

Gr=^S2*(® M ) +^£r 0M +^£o ( 48 ) 
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where 


82 Ka 

(49) 

ft-A-ey- a * 

(50) 

g o =A-0^0 Hr +A o ^ 

(51) 


Parameter g 2 represents the mechanical contribution to the strain energy release rate, 
while parameter g x is the coefficient of the cross term representing the coupling of 
hygrothermal and mechanical contributions, and g 0 represents the hygrothermal 
contribution. Equation 47 indicates that there could be a non zero strain energy release 
rate due to hygrothermal effects alone. 

Mode III Lav-up Identification 

The presence of extension-twist and bending-twist couplings in a laminate with 
arbitrary lay-up may induce mode I and mode II contributions to the total stain energy 
release rate in addition to the mode III contribution. The objective here is to identify the 
mode I and II contribution to the total strain energy release rate in order to minimize 
them. The actual evaluation of mode I and II components is not necessary. Instead, the 
closed form solution given in Equation 36 may provide a convenient alternative to 
identify the presence of mode I and II contributions. Substituting Equation 40 into 
Equation 49 yields 

82 = ~V 3 3 ,« - {*JV 31 + tfPjU ) - {<l>2,aVn + ^ 32 . a ) (52) 

The second and third terms in Equation 52 indicate the presence of mode I and II 
contributions due to extension-twist and bending-twist couplings, respectively. Both 
these terms should be zero for pure mode III lay-ups. If for a particular lay-up these two 
terms are very small compared to the first term, this lay-up may be considered as a 
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desirable mode III lay-up. The [90/(±45) n /(+45) n /90] s ECT lay-ups belong to this 
category. Two indicators are defined to help the identification of desirable mode ID lay- 
ups as follows 


_ = + (5 , 

^*33, a + + ( 02, 0^32 + 02 > ^32,o) 

_ ( l h ^_ 32 + 02^ 32,a ) 

“ 2 ~ ^3,a+{^n+^l,a) + {<l>2.a^2+02^2, a ) 

where and E 2 are extension- twist and bending- twist coupling indicators, 
respectively. Note that these indicators are percentages of the second and third terms to 
the total term in Equation 52. 


Results 

First, the present analytical solution was compared with Ref. 10 and 1 1 for the 
ECT lay-ups [90/(±45) n /(+45) n /90] s . The lamina properties are given in Table 1. The 
hygrothermal bending curvature and twist defined in Equation 33 vanish for these ECT 
lay-ups because of the laminate and sublaminate symmetry. For these ECT lay-ups, the 
present solution should be identical to that of Ref. 10. Tables 2 and 3 show the 
comparisons of torsional stiffness and strain energy release rate parameter (g 2 ) predicted 
by the present solution and Ref. 10 and 11. The present solution is identical to Ref. 10 as 
seen in Tables 2 and 3. Also observed is the accuracy of the simplified solution given in 
Ref. 1 1 for the ECT lay-ups. 

Altere d E CT La y -ups 

Second, the present analytical solution was used to investigate the residual 
thermal effect for several lay-ups selected by altering the [90/(±45) 3 /(+45) 3 /90] s ECT 
lay-up and are presented in Table 4. The lamina properties from Table 1 were used to 
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generate the analytical solutions and are presented in Table 5 for a delamination length 
a/b=03. The temperature change Ar=-156°C was used. The moisture weight gain AH 
was assumed to be zero to isolate the residual thermal effects. 

All these altered lay-ups have higher torsional stiffness (A,) than the original ECT 
lay-up when one compares the stiffness values in Table 5 to the corresponding value in 
Table 2. The thermal-mechanical coupling strain energy release rate parameter (gj) is 
zero, and the thermal parameter (go) is negligible compared to the magnitude (-1000 
N/m) of the mode III toughness obtained from experiments [9]. The small amount of g 0 is 
due to the asymmetry in the sublaminates. The symmetry of the entire laminate however 
results in the coupling parameter g l to diminish. The extension-twist coupling indicator 
(5j) vanishes for all the lay-ups because there is no extension-twist coupling in these lay 
ups. 


The bending-twist coupling indicator (E 2 ) is very small for these lay-ups as shown 
in Table 5 (0.6%<). Both lay-ups LI and L2 have six pairs of ±45° layers in each 
sublaminate, the bending-twist coupling indicator of LI is one order of magnitude 
smaller than that of L2 due to self-symmetry arrangemet of the ±45° layers in LI. The 
bending-twist coupling indicator is higher in lay-up L3 than that in lay-up L4 where the 
out-most 90° layer from the ECT lay-up has been changed to 45° and -45°, respectively. 
There are three and two pairs of ±45° layers in each sublaminate of lay-up L5 and L6, 
respectively. The bending-twist coupling indicator is higher in lay-up L6 than that in 
lay-up L5, but still less than 0.6%. All these lay-ups may be considered as desirable mode 
III lay-ups from the dominance of mode III strain energy release rate point of view. 
However, these lay-ups may suffer from distortion after curing process due to asymmetry 
in the sublaminates. 
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Antisymmetric Lav-ups 

The lay-ups studied above are all symmetric laminates and show little thermal 
effect on the strain energy release rate. The analytical model developed in Section 2 is 
not limited to symmetric laminates. Next, a class of asymmetric laminates are considered 
to investigate the residual thermal and moisture effects. These lay-ups are a class of 
antisymmetric laminates with [±(0/(0-9O) 2 /0)] n lay-ups. The sublaminates may be either 
symmetric (n=odd number) or antisymmetric (n=even number). This class of laminates 
exhibits extension-twist coupling, and does not distort after curing. 

Results for lay-ups 0=30° and n=l to 4 are presented in Table 6 for a normalized 
delamination length a/b=03. The moisture weight gain, A H, was assumed to be zero 
when the results in Table 6 were generated. Both torsional stiffness, strain energy release 
rate parameters and extension-twist coupling indicator are presented in the table. The 
bending-twist coupling indicator vanishes because there is no bending-twist coupling. 
The thermal parameter go is zero for all the lay-ups due to the symmetry or antisymmetry 
of the sublaminates. The presence of the coupling parameter g! and the extension-twist 
coupling indicator are strong when the sublaminates are symmetric (n=l and 3). 
However, the significance of parameter g t and mode I indicator reduces as n increases. 
The moisture effect may also be strong in these lay-ups. Fig. 4 shows the moisture effects 
on the coupling parameter g x as a function of delamination length a/b for n=3 with a 
fixed temperature change of A7=-156 °C. As the percentage moisture weight gain AH 
increases from zero, parameter gj reduces and becomes zero (canceling the thermal 
effect) when A//=0.0208 irrespective of the delamination length. Further increase of AH 
changes the sign of g t . 

When the sublaminates are antisymmetric (n=2 and 4), both the thermal- 
mechanical coupling parameter g x and the extension- twist coupling indicator are very 
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small. From this point of view, these lay-ups could be candidate mode III lay-ups to 
characterize the mode III fracture toughness between two layers with different fiber 
orientations. Additional results for 0=20°, 40°, 45°, 50°, 60° and 70° when n=4 are 
presented in Table 7 along with 0=30°. The torsional stiffness X, strain energy release 
rate parameter g 2 and the extension-twist coupling indicator are symmetric with 
respective to 0=45° as seen from Table 7. Both parameter g t and the extension-twist 
coupling indicator are very small and become identical zero when 0=45°. 

The torsional stiffness (X) and strain energy release rate parameter (g 2 ) as a 
function of normalized delamination length (a/b) are plotted in Fig. 5 and 6 for 0=30°, 
45°. The torsional stiffness decreases as the delamination length increases, and reaches 
the torsional stiffness of the two sublaminates as the delamination completely separates 
the laminate in to two pieces. The middle portion of the curve is almost linear. The curve 
becomes flat near both ends (a/Z>=0 and 1). This indicates that both very small and very 
large delamination lengths are not desirable for mode III delamination growth because 
the available strain energy release rate to drive the delamination is proportional to the 
derivative of the torsional stiffness. Fig. 6 is essentially the negative value of the 
derivative of the torsional stiffness. It is seen from Fig. 4 through 6 that the analytical 
solution gives meaningful results for the entire range of normalized delamination length 
(a/b=0 to 1). 

Concluding Remarks 

A shear deformation theory including hygrothermal effects is developed for the 
analysis of either symmetric or asymmetric laminates with mid-plane edge delamination 
under torsion loading. The torsional stiffness and total strain energy release rate are 
obtained in closed form. The hygrothermal effects on the strain energy release rate are 
identified by hygrothermal-mechanical coupling and pure hygrothermal contributions. 
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Two indicators are defined to identify the presence of mode I and mode II contributions 
in the total strain energy release rate. 

The analysis of the [90/(±45) n /(+45) n /90] s ECT mode III lay-ups indicates that 
there are no hygrothermal effects on the strain energy release rate because of the 
laminate and sublaminate symmetry. Identical results to those obtained in Ref. 10 are 
observed for these ECT lay-ups. A further study of some altered ECT lay-ups reveals 
that these lay-ups have negligible residual thermal and bending-twist coupling effects 
even when the sublaminates above and below the delamination are not symmetric about 
their own mid-planes. However, these lay-ups may suffer from distortion after the curing 
process. 

Another set of lay-ups investigated is a class of antisymmetric laminates with 
[±(0/(0-9O) 2 /0)] n lay-ups There is no bending-twist coupling effect in these lay-ups. 
When n takes on odd numbers (1 and 3), both the thermal and extension-twist coupling 
effects on the total strain energy release rate may be strong. The relative significance of 
the thermal and extension- twist coupling effects reduces as n increases. The moisture 
effect may also be strong in these lay-ups. However, the moisture effect tends to reduce 
the thermal effect, and may cancel the thermal effect at some value of the percentage 
moisture weight gain irrespective of the delamination length. 

When n takes on even numbers (2 and 4), both the hygrothermal and extension- 
twist coupling effects may be neglected. From this point of view, these lay-ups are 
desirable mode III lay-ups to characterize the mode III fracture toughness between two 
±0 layers. In particular, when 0=45°, both the hygrothermal and extension-twist 
coupling effects vanish irrespective of n. 
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Appendix : Definition of Parameters 

In this appendix, all the intermediate parameters needed to evaluate the torsional 

stiffness and strain energy release rate are summarized according to their sequence of 

appearance in the text. 

The matrix [R j ] and vectors {fl j } and {fl j } in Equation 22 are given by 
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The integration constant, P , in Equation 26 is given by 
/j=/ jk tf, j,k =1,2,3 

The explicit expressions of /j k are given in the following 
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The matrix ['Ey] and vector { 4 / | rl } in Equation 32 are given by 
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TABLE 1-- Elastic properties and configurational parameters of Carbon/Epoxy composite 

[9] 


Properties* 

Configurational Parameters 

E n =165 GPa 

<Xi=-0.3 x 10' 6 /°C 

b= 38.1 mm (Width) 

E 22 =E 33 =10.3 GPa 

a 2 =30 x 10- 6 /°C 

r=0.13 mm (Ply thickness) 

GPa 

Pi=0 

/i=Half laminate thickness 

// 23 =5.5 GPa 

P 2 =0.2 


Vi2 = Vi3=0.28 

Ar=-156°C 


v 23 =0.28 




* These thermal and moisture expansion coefficients are adopted from a similar 
carbon/epoxy (IM6/SC1081) material [15] because such information was not available in 
Ref. 9. 


TABLE 2-- Comparison of torsional stiffness and strain energy release rate parameter g 2 
for the ECT lay-up [90/(±45) n /(+45) n /90] s with n=3 


n=3 

X (Nm 2 ) 

g, (Nm) 

a/b 

Present 

Ref. 10 

Ref. 11 

Present 

Ref. 10 

Ref. 11 

0.1 

17.49 

17.49 

17.53 

426.4 

426.4 

424.4 

0.3 

14.06 

14.06 

14.11 

451.0 

451.1 

449.7 

0.5 

10.60 

10.60 

10.66 

450.1 

450.2 

448.8 

0.7 

7.244 

7.244 

7.311 

434.6 

434.7 

433.3 

0.9 

4.434 

4.434 

4.510 

223.7 

223.9 

223.1 
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TABLE 3-- Comparison of torsional stiffness strain energy release rate parameter g 2 for 
the ECT lay-up [90/(±45) n /(T45) n /90] s with n=4 


n=4 

X (Nm 2 ) 

g 2 (Nm) 

a/b 

Present 

Ref. 10 

Ref. 11 

Present 

Ref. 10 

Ref. 11 

0.1 

37.19 

37.19 

37.24 

859.7 

859.8 

855.9 

0.3 

29.82 

29.82 

29.88 

985.9 

986.0 

985.7 

0.5 

22.27 

22.27 

22.36 

976.5 

976.6 

976.3 

0.7 

15.13 

15.13 

15.23 

890.7 

890.8 

890.5 

0.9 

9.852 

9.852 

9.964 

360.0 

360.2 

359.6 


TABLE 4— Altered ECT lay-ups and their representative codes. 


Code 

Lay-up 

LI 

[(±45) 3 /(+45) 3 /0/90] s 

L2 

[(±45)6/ 0/90] s 

L3 

[45/(±45) 3 /(+45) 3 /90] s 

L4 

[-45/(±45) 3 /(+45) 3 /90] s 

L5 

[(±45 2 ) 3 /0/90] s 

L6 

K±45 3 yo/90], 
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TABLE 5-- Torsional stiffness, total energy release rate parameters and bending-twist 
coupling indicator for the altered ECT lay-ups at a/b=03 


Code 

A.(Nm2) 

g? (Nm) 

gi (N) 

go (N/m) 

s 2 % 

LI 

16.39 

567.4 

0 

0.040 

-0.002 

L2 

16.31 

561.7 

0 

0.040 

-0.057 

L3 

16.06 

526.5 

0 

0.278 

-0.323 

L4 

16.54 

531.0 

0 

0.277 

-0.185 

L5 

15.91 

549.5 

0 

0.040 

-0.234 

L6 

15.24 

529.2 

0 

0.041 

-0.544 


TABLE 6- Torsional stiffness, total energy release rate parameters, and extension-twist 
coupling indicator for [±(30/-60 2 /30)] n lay-ups at a/b=Q3 


n 

X (Nm 2 ) 

g 2 (Nm) 

g. (N) 

go (N/m) 

E 1 % 

1* 

.247 

5.244 

8.078 

0 

-64.5 

3 

8.49 

247.7 

33.05 

0 

4.1 

2* * 

2.39 

74.98 

-0.5 x 10- 5 

0 

0.3 x 10- 5 

4 

19.5 

599.4 

-0.042 

0 

0.003 


* [(30/-60 2 /30) / (-30/60 2 /-30)] 

T 

**[(30/-60 2 /30/-30/60 2 /-30)/(30/60 2 /30/-30/60 2 /-30)] 

Note that T indicates the interface with free edge delamination. 
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TABLE 7- Torsional stiffness, total energy release rate parameters, and extension-twist 
coupling indicator for [±(0/(9O-0) 2 /0)] 4 lay-ups at a/b=03 


0 

A,(Nm 2 ) 

g 2 (Nm) 

g,(N) 

5! % 

20 

12.6 

347.9 

-0.004 

0.0003 

30 

19.5 

599.4 

-0.042 

0.003 

40 

24.3 

745.1 

-0.046 

0.002 

45 

25.0 

765.0 

0 

0 

50 

24.3 

745.1 

0.046 

0.002 

60 

19.5 

599.4 

0.042 

0.003 

70 

12.6 

374.9 

0.004 

0.0003 
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FIG. 1- A laminate with mid-plane edge delamination under torsion loading. 
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FIG. 2— A typical x= constant plane of the delaminated laminate. 
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FIG. 3— Notation and sign convention for stress resultants and moment couples. 
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FIG. 4— Coupling strain energy release rate parameter, g lf as a function of normalized 
delamination length, a/b , for [±(30/-60 2 /30] 3 lay-up. 
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